DELPHI Collaboration 



DELPHI 98-42 PHYS 772 
27 April, 1988 



Study of b — > uiv Decays 
with an Inclusive Generator 



M. Battaglia 

Dept. of Physics, 
University of Helsinki (Finland) 



Abstract 

A generator for inclusive b — ► uiv decays has been developed. Different prescriptions 
have been applied to describe the kinematics of the b quark inside the hadron and 
the generation of the hadronic final states. Results are presented with particular 
attention to the invariant mass of the hadronic system recoiling against the lepton 
and its resonance decomposition. These studies are of special relevance for the 
extraction of \V u b\ from semileptonic B decays at LEP and at B factories. 
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1 Introduction 



The measurement of the branching ratio for the decay b — > ulv provides the cleanest way 
to determine the \ V u b\ element in the CKM mixing matrix. Evidence of the non-zero value 
of \V u b\ has been first obtained by both ARGUS and CLEO by observing leptons produced 
in B decays with momentum exceeding the kinematical limit for b — > civ transitions ]l|, . 
The extraction of \V u b\ from the yield of leptons above the b — > civ endpoint is subject 
to large systematical uncertainties. More recently, exclusive B — > irlv and B — > plv 
decays have been observed by CLEO and their rates measured || |J. Still the derivation 
of \V u b\ from exclusive semileptonic decays, contrary to the case for \V c b\ in B — > D*lv, 
has significant model dependence. 

The extraction of \V u b\ from the shape of the invariant mass of the hadronic system 
recoiling against the lepton in b — > ulv transitions was proposed several years ago || and 
it has been recently the subject of new studies || [7J. The proposed method starts from 
the observation that the hadronic system recoiling against the lepton in the decay has 
invariant mass lower than the charm mass for the majority of b — > ulv decays. The model 
dependence in predicting the shape of this invariant mass distribution has been claimed 
to be under control within about 10 — 15 % if b — > ulv decays can be distinguished from 
the b — > civ ones for masses of the hadronic system up to cut values close enough to the 
D mass 0. This corresponds to a model uncertainty of 5-7% on the extraction of \V u b\. 

The predicted shape of the invariant mass distribution depends mainly on the kine- 
matics of the heavy and spectator quarks inside the B hadron and on the quark masses. 
However from the experimental point of view, the hadronisation process, transforming 
the uq system into the observable hadronic final state, represents a significant source of 
additional model uncertainties. 

Several models have been proposed to describe both steps of the decay process. This 
note summarizes the results of a study performed by developing a dedicated b — > ulv 
decay generator, BTOOL. This generator implements different prescriptions for the initial 
state kinematics and the resonance decomposition of the hadronic final states. Its results 
are used to define the model dependent systematics in the extraction of \V u b\ from the 
hadronic mass spectrum in semileptonic B decays. 

2 The BTOOL Generator 

The decay generator provides the four momenta of the stable decay products in b — > ulv 
transitions. This requires a model for the kinematics of the b and spectator q quarks inside 
the B hadron, the description of the Q 2 distribution of the virtual W and of the kinematics 
in the b — > uW and W —>■ Iv decays, and finally a model for the hadronisation of the uq 
system. In the following the implementation of the decay in the generator is discussed. 
In obtain to study model dependences and systematic effects, different prescriptions have 
been adopted. 

2.1 ACCMM Model 

In the ACCMM model || the B hadron consists of the b quark and the spectator q quark 
moving back-to-back in the B rest frame. Their momenta p are distributed according to 
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the gaussian distribution: 

where the width pp is known as Fermi motion and represents a parameter in the model. 
The normalisation is chosen such that 

/ dpp 2 <P(p) = l. (2) 
Jo 

The choice of the pf parameter and of the mass of the spectator quark m q are discussed 
in details in the next section. 



2.2 Parton Model 

An alternative picture of the b quark kinematics has been proposed as an application of 
the parton model to heavy quark decays |J. In this model the decay is considered in 
a frame where the B hadron moves with large momentum (infinite momentum frame or 
Breit frame). In this frame the b quark behaves as a free particle earring a fraction z of 
the B momentum, pt, = zps- The functional form of f(z) can be extracted from the b 
fragmentation function since the probability of finding the b quark carrying a fraction z of 
the hadron momentum corresponds to the probability for a b quark to produce a B hadron 
with a fraction z of its energy. This is usually described by the Peterson fragmentation 
function |10| that has the form: 



Hz) - NZ{1 ~ Z? (3) 

where e& is a free parameter. In this way, the parton model offers an advantage since 
the kinematics of the b quark is described by a function that can be directly related with 
experimental data on b fragmentation. 

2.3 QCD Universal Structure Function 

Recently there has been progress in defining the Fermi motion in the framework of 
QCD JTT], |TJ, [TBJ. This has been achieved in terms of a universal structure function 



describing the distribution of the light-cone residual momentum of the heavy quark inside 
the hadron. At leading order and in the large m^ limit, the light-cone residual momentum 
k + can be expressed as the difference between the b quark pole mass and its effective mass 
ml inside the hadron: ml = m^ + k + . As a consequence k + < A = m# — m^. An ansatz 



for the shape of the universal structure function has been suggested [13] in the form 



f{z) = z\l-cz)e- cz (4) 

where z = 1 — and the coefficients a and c depend on the values of A and of the kinetic 
energy operator as discussed in the next section. 
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2.4 Decay kinematics 

The b quark decays as b — > Wu, the W and w being emitted back to back in the b rest 
frame with an isotropic distribution of their emission angle w.r.t the b direction. The 
virtual W is characterized by an effective mass Q 2 . The kinematics of the b and W 
decays correspond to that for two body decays. Therefore the model dependence that 
propagates to the final state kinematics depends on the choice of the values of pp and m q 
and of the Q 2 distribution. By defining x 2 = Q 2 /m 2 the differential decay rate can be 
expressed as: 



dx 



G 2 F ml\V uh \ 

1927T 3 



2a 

<Fo{x) - -^F x {x)) 



(5) 



The functions F (x) and Fi(x) describe the tree- level contribution and the QCD cor- 



rections terms. Following Ref. |TJ[, they can be written, in the limit m u 



0, as: 



F Q (x) = 2(1 



x 2 ) 2 (l + 2x 2 ) 



and 



= (tt 2 + 2S 1A (x 2 ) - 2S ltl (l - x 2 )) + 8x 2 (l - x 2 - 2x A )ln(x) 
+2(1 - x 2 ) 2 (5 + 4x 2 )/n(l - x 2 ) - (1 - x 2 )(5 + 9x 2 - 6x 4 ) 



(6) 



(7) 



where Si t i(x) is the Nielsen polylogarithm. The resulting Q 2 distribution is shown in 
Figure 1. 




22.5 25 

q 2 (QeV 2 ) 



Figure 1: The q 2 distribution in the decay b — > utv generated according to Eq. (5). 



The virtual W is forced to decay to the lepton-neutrino pair. Due to the spin of the 
W, its decay is not isotropic and the angular distribution of the lepton varies as 1 + cos 2 9. 
In computing the lepton and neutrino momenta in the W rest frame the lepton mass is 
taken into account. After the decay the lepton and neutrino are boosted back to the B 
rest frame. The u quark is also boosted to the same system, and at this point the hadronic 
system is treated. 
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2.5 The hadronic system 



The energy and invariant mass of the hadronic system correspond to those of the uq 
quark pair. The energy and invariant mass at the parton level can be compared with 
predictions obtained in QCD and heavy quark expansion. The observable final states 
are then generated by applying different prescriptions for describing the evolution of the 
quark fragmentation as discussed in the next section. These include a fully inclusive 
hadronisation scheme according to the JETSET parton shower model |15| and exclusive 
descriptions of the final states in the b — > ulv transition. 



2.6 The JETSET 7.4 interface 

The generator can be used to produce individual B meson decays for dedicated studies. 
Timings for the generation of individual decays on different platforms are given in Table 1. 
BTOOL is also interfaced as an option of the LUDECY subroutine of JETSET to handle 
b — > ulv decays in the generation of e + e~ — > Z°/7 — * bb events. The BTOOL generator 
can be activated by setting the MDME(IDC,2) flag in the LUDAT3 common block of 
JETSET, where IDC refers to the b — > ulv decay channel. 



Table 1 : Timing of the decay generator for different platforms 



Platform 


OS 


secs./lOOO decays 


HP 9000/778 
Pentium- 133 


HP-UX 10.20 
Linux 2.0.27 


0.52 
1.85 



3 Results and Comparisons 

The focal interest in the simulation study of b — > ulv decays is to determine the invariant 
mass spectrum of the hadronic system and its correlation with the lepton energy and then 
to define the uncertainties of these distributions due to the choice of the model and of 
its input parameters. Three different models have been applied for the definition of the 
kinematics of the b and spectator q inside the hadrons and two models for the generation 
of the hadronic final states. The criteria chosen for the input parameters and their range 
of variation are discussed in the following subsection. The results for the hadronic system 
mass and multiplicities are presented in 3.2 and 3.3. 

3.1 The Choice of Parameters 

The ACCMM model introduces two free parameters: i) the Fermi momentum pp and 
%%) the spectator quark mass m sp . There have been attempts to extract the values of 
these parameters from fits to experimental observables such as the momentum spectrum 
of leptons from b — ► civ decay |]I6| [lTj , the photon spectrum in b —>■ s'j decay |18[] and the 
J/ip momentum distribution in B —>■ J/ip X [Tj|. Most of these determinations point to 
a value of p F ~ 0.5 GeV/c (Table 2). 
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Table 2: Experimental estimates of the pp value in the ACCMM model 



Channel 


p F (GeV/c) 


m sp (GeV/c 2 ) 


Ref. 


b — > civ 
b — > civ 
b — > 57 


0.27 ± 0.04 

u - 0i -0.07 

0.45 
0.57 


0.30 
0.0 
0.0 
0.15 


[ 

1 


16 
17 
18 
19 





However it has been pointed out that the value of pp obtained in a fit to b — > c 
transitions may be not appropriate in the description of b — > w decays |p0| . At the 



same time it has also been shown that the ACCMM model is consistent with the QCD 
description of b —>■ ulv and b —>■ s'j transitions and that the corresponding parameters 
can therefore be related. This is the direction followed in this study. 
The effective b quark mass m& depends on pp and 



™Z = ml(p b ) = m| + m 2 sp - 2m B ^p 2 h + m 2 sp (9) 

where pb is the momentum of the heavy quark in the hadron and is the B hadron 
mass. The value ofmp can be taken as a parameter, tunable such that the ACCMM model 
corresponding to an average b mass < m& > can be compared with theory predictions 
obtained for a given value of m b . Estimated values for the b quark pole mass are in the 
range 4.72 GeV/c 2 < m b < 4.92 GeV/c 2 Q. 



The value of pp is proportional to the average kinetic energy of the b quark in the 
hadron since: 

r+oo 3 

< vl >= / dp b pi (<j){p b ) pi) = -p 2 F (10) 
Jo 2 

where <fi(p) is given by Eq. (1). Through the two above equations the ACCMM model 
parameters are related to those of the QCD description of the heavy quark inside the 
hadron. In this framework < p\ > corresponds to the value of the expectation value [x\ 
of the kinetic operator. The value pp = 0.5 GeV/c corresponds to < pi > = 0.37 GeV 2 . 
Estimates of < p\ > have been obtained both from theory and fits to measured spectra 
in B decays as discussed below. 

For the Parton Model, the fragmentation function for b quarks has been measured at 
LEP. Averaging over the ALEPH, DELPHI and OPAL results, the fraction of the b quark 



energy taken by the beauty hadron is < ig >= 0.702 ± 0.008 [21|. Also the observed 



shape in the preliminary DELPHI analysis |22[] was compatible with that of the Peterson 
function. These results point to a value for the e b parameter of e b = 0.0040. In the 
simulation of decays with the parton model, the spectator quark mass mass m q was set 
to zero and the b quark mass was varied in the range 4.72 GeV/c 2 < mj < 4.92 GeV/c 2 
as for the ACCMM model. It is interesting to point out that, by using the central value 
for e b , the parton model gave < pi >= 0.35 GeV 2 for m b = 4.72 GeV/c 2 , consistent with 
the value obtained in the ACCMM model for pp = 0.5 GeV/c. 

The use of the QCD universal structure function f{k + ) allows a consistent comparison 
with the results obtained in the framework of QCD and Heavy Quark expansion. The 



5 



normalised moments a r 



j* of f(k+), given by 



111 



relate the function parameters with that of the theory. In particular the first two moments 
define the function normalisation and the third is proportional to the expectation value 
of the kinetic energy operator | T3fl : 



(1-2 



a = 1 

a± = 

3/4 
A 2 



(12) 
(13) 

(14) 



These relationships define the values of the parameters a and c in Eq. 4 as a function 
of the values of /i 2 and A. There have been several evaluations of /z 2 and a selection of 
recent results is given in Table 3. Results are scheme dependent and, depending on the 
method used in their derivation, they point to the values of of 0.4 GeV 2 or 0.2 GeV 2 . 



Table 3: Estimates of /x 2 



Method 


A (GeV 2 ) 


Ref. 


QCD Sum Rules 


0.50 




23 




Variational Method 


0.44 




24 




Relativ. Potential 


0.46 




25 




Virial Theorem 


0.40 - 0.58 


[ 


26 




Virial Theorem 


0.15 ± 0.03 


|27fl 


civ 


0.19 ± 0.10 




28 


] 


b^civ 


0.14 ± 0.03 




29 




b — > 57 


71 + 116 
u - ' i -0.70 




30 





The two values of /i 2 = 0.2 GeV 2 and 0.4 GeV 2 were chosen while A ranged be- 
tween 0.36 GeV/c 2 and 0.56 GeV/c 2 for 4.72 GeV/c 2 < m h < 4.92 GeV/c 2 and m B = 
5.28 GeV/c 2 . The masses of the light quarks, u and spectator quark, were set to zero. 

3.2 The Hadronic System Mass 

The precise determination of the fraction of b —>■ uiv transitions yielding an hadronic 
system with mass Mx below a given M cut value is crucial in the estimation of \V u b\. As 
discussed above, at the quark level the hadronic mass corresponds to the mass of the 
uq system. This depends mainly on the masses and motion of the heavy and spectator 
quarks in the B hadron. On the contrary the experimentally measureable hadronic mass 
is strongly affected by the resonant decomposition of the hadronic final states. In the 
simulation the hadronic system is analyzed in two steps. 

Firstly, the uq pair is analyzed. This gives the hadronic mass distribution before 
resonant states are taken into account (see Figure 2) and can be compared with that 
computed using QCD and heavy quark expansion 0. 
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u q Mass (GeV/c 2 ) u q Mass (GeV/c 2 ) 



Figure 2: The invariant mass distribution of the uq system obtained using the ACCMM 
model (left) and the QCD universal function (right) corresponding to nib = 4.72 GeV/c 2 
and <pl> = OA GeV 2 . 

Results are expressed in terms of the fraction F u (M cut ) of the b — > utv transitions 
resulting in a mass of the uq pair below a given cut value. The kinematics of the b and 
spectator quark have been defined using the ACCMM model, the QCD universal structure 
function and the Parton Model. In order to compare the results, parameters have been 
chosen such that they correspond to nib = 4.80 ± 0.10 GeV/c 2 and 0.2 GeV 2 < < 
pi > < 0.4 GeV 2 as discussed above. 



Table 4: Fraction F u (M cut ) ofb^ utv decays with M x < M cut for < p 2 b > = 0.4 GeV/c 2 . 



ni h 


M cut 


F u (M cut ) 


F u (M cut ) 


F u (M cut ) 


F u (M cut ) 


(GeV/c 2 ) 


(GeV/c 2 ) 


ACCMM 


QCD Funct. 


Parton 


Ref. § 


4.92 


1.25 


0.48 


0.50 


0.40 


0.45 




1.50 


0.65 


0.60 


0.62 


0.60 




1.75 


0.76 


0.80 


0.84 


0.78 


4.82 


1.25 


0.50 


0.57 


0.45 


0.55 




1.50 


0.65 


0.72 


0.74 


0.72 




1.75 


0.78 


0.83 


0.89 


0.85 


4.72 


1.25 


0.54 


0.65 


0.47 


0.68 




1.50 


0.68 


0.75 


0.71 


0.81 




1.75 


0.80 


0.85 


0.87 


0.89 



Table 4 summarizes the results for different choices of the input parameters. The 
first observation is that the universal function implemented in BTOOL reproduces the 
prediction from QCD and heavy quark expansion from 0. Further it has been found 
that also the ACCMM model reproduces these results to better than 15%, for equivalent 
values of m& and < p 2 >. This study confirms that the sensitivity to the value of the b 
quark mass is significant in the region of low hadronic invariant masses. In this region also 
the model dependence is more pronounced. In order to estimate the overall uncertainty 
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in the estimate of the fraction of b — > ulv decays with hadronic mass below a given cut 
value M cut , the different sources of systematics have been combined. The chosen range of 
variation of the parameters is ±a(mb) = ± 0.10 GeV/c 2 and ±<r(< pi >) = ± 0.1 GeV 2 . 
The corresponding realtive systematic errors are summarized in Table 5. 



Table 5: Estimate of the relative systematic uncertainties on F u (M cut ) 



Source 


M cut = 1.25 


M cut = 1.50 


M cut = 1.75 


m b 


0.13 


0.10 


0.05 


<Pb > 


0.08 


0.04 


0.02 


Model 


0.08 


0.04 


0.04 


Total 


0.17 


0.12 


0.07 



Secondly, the hadronic final states corresponding to a given mass and energy of the 
uq pair can be predicted by a variety of methods ranging from the fully inclusive quark 
fragmentation approach to exclusive models. 
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Figure 3: The distribution of the u quark energy versus the lepton energy (a) and the 
physical region in the q 2 versus 2m#go plane (b) for inclusive b — » utv decay. 

At large enough recoil u quark energies, the uq system moves away fast and this 
picture is similar to that of the evolution of a jet initiated by a light quark q in e + e~ — > qq 
annihilation. This is simulated by first arranging the uq system in a string configuration 
and then making it fragment according to the parton shower model. Exclusive models 
compute the decay amplitudes from the heavy-to-light form factors and the quark hadronic 



wave functions. The so-called ISGW2 model |31|j approximates the inclusive b — > ulv 
decay width by the sum over resonant final states, taking into account leading corrections 
to the heavy quark symmetry limit. 

The range of applicability of the inclusive and exclusive models is restricted to par- 
ticular regions of the accessible kinematic configurations. Figure 3 a) shows the physical 
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region in the q 2 versus 2 tub qo plane where q 2 is the effective mass of the virtual W and 
go its energy in the B rest frame. In this plot states of equal hadronic invariant mass Mx 
correspond to lines M\ = m? B — 2rriBqo + q 2 ■ Systems of large invariant mass correspond 
to low q 2 and go values, i.e. large u recoil as shown in Figure 3 b). In these cases the 
u quark energy is typically large enough compared with that of the spectator quark that 
the analogy with jet fragmentation is justifiable. Conversely at low u recoil energy, i.e. 
close to the upper kinematical limit in Figure 3 a), the relative momentum of the uq pair 
is small and they are therefore likely to form a bound state. 




Figure 4: Comparison of the invariant mass of the hadronic final state from the ISGW2 
model (upper left) and from parton shower fragmentation (upper right). In the lower plot 
the corresponding fractions F u (M cut ) of the b — > uiv transitions resulting in a mass below 
a given cut value M cut are shown as a function of M cut . The corresponding predictions 
for the invariant mass of the uq pair are shown with the continuous line (jjffl) and the 
dash histogram (BTOOL). The CLEO results for the branching ratios for the n£i? and 
piv channels are also shown for comparison. 

A satisfactory description of the inclusive b — *> uiv decay must take these character- 



istics into account. Hybrid models have been proposed for this purpose [|3^]. The main 
feature of an hybrid model is to define a kinematical region in which the exclusive model 
is valid and its complement that can be treated by an inclusive fragmentation model. The 
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two regions must be chosen in order to have well behaved matching conditions for a set 
of relevant kinematical variables. Further constraints can be derived from the branching 
ratios for B — > tt£u and B — > ptv decays measured by CLEO [[3], |J. Since the aim of this 
study is the definition of the systematic uncertainties in the description of the hadronic 
mass spectrum, results have been derived for the two extreme cases of fully inclusive and 
exclusive models. In the inclusive model, the probabilities for generating light vector and 
axial resonances have been tuned in JETSET in order agree with the measured rates in 
Z° decays. For the exclusive model the ISGW2 model has been used. The results are 
presented in Figure 4 in terms of the fraction of b — > ulv decays with the hadronic final 
state below a given mass value. The comparison of the predictions from the two models 
shows significant differences in the predicted mass spectra due to the relative importance 
of resonant and non-resonant final states. For the region of M cut > 1.5 GeV/c 2 , the rela- 
tive difference of the two models correponds to 10 - 15% showing that the hadronic system 
fragmentation introduces an uncertainty of the same order as that from the 6-quark mass 
and the heavy hadron kinematics. 

3.3 The Hadronic System Multiplicity 

An additional source of uncertainty in the modelling of the decay arises from the multi- 
plicity of the hadronic system. This is of special relevance since the efficiency for reson- 
structing the decay depends on this multiplicity. In order to study this uncertainty, the 
different prescriptions for describing the final states discussed above have been analyzed 
in terms of the resulting decay multiplicity. 




01 2345678 



E cms (GeV) 



Figure 5: Comparison of the generator predictions with experimental data on charged 
multiplicity in low-energy e + e~ collisions. The points represent the experimental data, 
the long thick line the best fit from Ref. |33|/ and the shorter thinner lines the fits to the 
generator results for values ofpF ranging from 0.0 GeV/c to 0.35 GeV /c. 
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The charged multiplicity of the hadronic system from the b — > utv decay can be 
compared with one half of that of a qq event at y/s = 2Eh a d- The data on the event 
charged multiplicity from Adone and Mark II at 2 GeV < y/s < 8 GeV can be described 
with a function < n c h >= a + b In s with a = 2.67 ±0.04 and b = 0.48 ±0.02 [33| shown by 
the long thick line in Figure 5. These data have also been compared with the results of the 
simulation where a u quark of energy y/s/ 2 is paired with the spectator quark. The quark 
is either given Fermi motion according to the value of pp but no transverse momentum, or 
kept at rest (i.e. pp = 0). The results are shown by the shorter thinner lines in Figure 5. 
The parton shower model reproduces reasonably well both the multiplicity and its scaling 
with the quark energy. The best agreement with the data is obtained by imposing pp = 
0.2 GeV/c which gives a fit with a=2.71 and 6=0.44. 

Multiplicities in semileptonic B decays have also been predicted using the quark- gluon 
string model (QGSM) which also reproduces fairly precisely the same data on the charged 
event multiplicity in low energy e + e~ collisions |Q. The multiplicities obtained by the 
decay generator using the hybrid model are compared with the predictions from QGSM 
and IGSW2 in Table 6. The average charged multiplicity in the decay < n c h > agrees 
for the three models within ± 0.16. This multiplicity is also quite close to that measured 
for D meson decays |35[ , showing that the reconstruction of the hadronic system may be 
performed with comparable efficiency in semileptonic b —>■ u and b —>■ c decays. 



Table 6: Decay multiplicities from the event generator compared with QGSM and IGSW2 
models 



Final State 


Simulated 


QGSM @] 


ISGW2 





TT* 1 7T U 


.050 


.020 


.044 


7T ± 2 7T° 


.016 


.002 


.010 


7T ± 3 7T° 


.011 


.009 


.012 


7T ± 4 7T° 


.005 


.004 


.004 


1 7T ± 7T° 


.101 


.060 


.090 


1 7T ± 1 7T° 


.124 


.140 


.135 


1 TT* 2 7T° 


.037 


.064 


.016 


1 7T ± 3 7T° 


.008 


.025 


.008 


1 7T ± 4 7T° 


.003 


.003 




2 7T ± 7T° 


.117 


.080 


.090 


2 7T ± 1 7T° 


.176 


.175 


.167 


2 7T ± 2 7T° 


.038 


.095 


.048 


2 7T ± 3 7T° 


.004 


.015 


.005 


3 7T* 7T° 


.053 


.070 


.054 


3 7T ± 1 7T° 


.062 


.098 


.061 


3 7T ± 2 7T° 


.021 


.036 




4 7T ± 7T° 


.012 


.050 


.006 


4 7T ± 1 7T° 


.017 


.020 


.003 


< Weft > 


1.91 


2.06 


1.70 
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4 The Lepton Spectrum 



As already mentioned, the lepton spectrum is sensitive to the mass of the quark produced 
in the semileptonic b decay. While the lepton yield in the region of lepton energies 
above the kinematical limit ~ M f M ^ D for b — > civ transitions is subject to significant 
model dependences, a combined study of the mass of the hadronic system Mx and the 
energy of the lepton in the B rest frame Eg may allow an extraction of the \V u b\ with 
good sensitivity and improved control of the systematics. Figure 6 shows the correlation 
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Figure 6: Invariant mass of the hadronic system vs. lepton energy in b — > utv decays for 
the uq system (left) and the hadronic final state (right). 

between the values of Eg and Mx before and after the hadronisation of the uq system. 
By selecting decays with low hadronic invariant mass, the lepton spectrum is depleted in 
its lower end without significantly affecting the region of lepton energies above 1.5 GeV/c 
that are relevant for separating b — > uiv from b — > civ decays. 



5 Conclusion 

A generator for inclusive b — > uiv decays has been developed and used for studying 
the invariant mass and resonance decomposition of the hadronic system produced in the 
decay. These studies are of special relevance for the extraction of \V u b\ from semileptonic 
B decays at LEP and at the B factories. The low invariant mass of the hadronic system 
emitted in these decays can be used to separate them from the CKM favoured b —>■ civ 
transitions. Different models for defining the kinematics of the heavy quark inside the 
hadron and the uq system hadronisation have been compared. Systematic uncertainties 
in the fraction of decays giving an hadronic system with mass below a given cut arise from 
the value of the b quark pole mass, the momentum distribution of the heavy and spectator 
quark inside the hadron and modelling of the uq hadronisation. This analysis confirmed 
the observation that model dependencies and these systematic uncertainties from the b- 
quark mass and the heavy quark kinematics can be kept at the 10% level, or below, if 
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the study of b — > uiv is performed including decays with hadronic final state masses 
up to ~ 1.6 GeV/c 2 or above. A comparable uncertainty arises from the hadronisation 
model when comparing an inclusive to a fully exclusive model. The combined analysis of 
the hadronic mass and lepton spectrum spectrum may provide an optimal separation of 
b — > utv from b — > civ decays. 
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